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Abstract. Rigorous estimates for the total - (kinetic) energy plus pressure 
- flux in IR 3 are obtained from the three dimensional Navier-Stokes equations. 
The bounds are used to establish a condition - involving Taylor length scale 
and the size of the domain - sufficient for existence of the inertial range and 
the energy cascade in decaying turbulence (zero driving force, non-increasing 
global energy). Several manifestations of the locality of the flux under this 
condition are obtained. All the scales involved are actual physical scales in R 3 
and no regularity or homogeneity /scaling assumptions arc made. 



1. INTRODUCTION 

One of the main features of Kolmogorov's empirical turbulence theory [TT | [T2 l [T0] 
is existence of energy cascade over a wide range of length scales, called the inertial 
range, where the dissipation effects are dominated by the transport of energy from 
higher to lower scales. Energy cascades have been observed in physical experiments, 
but theoretical justification of this phenomenon using equations of fluid motion, and 
in particular, the Navier-Stokes equations (NSE), remains far from being settled. 
The technical complexity of the NSE makes it difficult to establish conditions un- 
der which such cascades can occur. A particular problem is the possible lack of 
regularity of the solutions to the NSE, and thus choosing the right setting becomes 
crucial. (For an overview of various mathematical models of turbulence and the 
theory of the NSE, see, e.g., [9j |6] and O SI EE] , respectively.) 

The first studies in this direction were made in [5], where infinite-time averages 
of the Leray-Hopf solutions in the Fourier setting were used to establish a sufficient 
condition for the energy cascade. This condition, involving Taylor length scale, 
provided an inspiration for the sufficient condition (|4.13|) obtained in section 4. In 
contrast to [8], our goal was to work in physical space, dealing with actual length 
scales in R 3 rather than the Fourier wave numbers. 

In studying a PDE model, a natural way of introducing a concept of scale is to 
measure oscillations, i.e., (distributional) derivatives of a quantity with respect to 
the scale. 

Considering an L\ oc function / on a ball of radius 2R, -B(xo, 2i?), the physical 
scale R is introduced via bounds on the distributional derivatives of / where a 
test function i/j is a refined - smooth, non-negative, equal to 1 on B(x.o,R) and 
featuring optimal bounds on the derivatives over the outer i?-layer - cut-off function 
on -B(xo, 2R). (Uniformity in all scales dictates linearity of the length of the outer 
layer in R; hence -B(xq, R + R).) 
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More explicitly, 

1(^7, VOI < I \f\\D a ^\ < (c(a)-^\f\,r {a) ) 

JB(x ,2i?) V ™ 1 7 

for some c(a) > and p(a) in (0, 1). 

(An attempt to introduce a concept of scale via characteristic functions in place 
of smooth cut-off functions would lead to infinite concentration - delta functions - 
invalidating much of the desired calculus.) 

This approach has a similar flavor as introducing the Fourier scale |£| via 

(in the Schwarz space, and then by duality in the space of tempered distributions) . 

Let xo be in B(0, Rq) (Rq being the integral scale, B(0, 2Rq) C £1 where f2 is the 
global spatial domain) and < R < i?o- Define local - per unit of mass - (kinetic) 
energy, e and enstrophy, E, at time t, associated with the ball i?(x , R) by 

!(*)= /jH^dx, 



e x ,M 
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Ex 



R (t) = J |V® u\ 2 cj)dx 



where 4> — V ip an d V an d ip are refined cut-off functions in time and space, respec- 
tively (for some | < 6 < 1). 

A total flux - (kinetic) energy plus pressure - through the boundary of a region 
D is given by 

/ (^' U ' 2 +P ) U ' nds = / [( u ' v ) u + v p] ' udx 

dD D 

where n is an outward normal. Considering the NSE localized to B(x ,2R) - and 
utilizing V • u = - leads to a localized flux, 

*xo,fl(*) = f Q|u| 2 +^u-V ( /)rfx = -y [(u-V)u + Vp] -u0dx. 

Since ^ can be constructed such that V<p = r]Vtl> is oriented along the radial 
directions of -B(xo, 2R) toward the center of the ball, ^xo.-R represents the flux into 
B(x ,R) through the layer between the spheres S(x ,2R) and S(x ,R) (V</> = 
onB(x ,i?)). 

A more dynamic physical significance of the sign of ^xoj-R can De seen from the 
equations: multiplying the NSE by tpu and integrating over B(x ,2R) (formally, 
assuming smoothness) leads to 

(1.1) j t J i|u| 2 i/.dx = $ X0:K + i/ J Au-U7/»dx. 

Plainly, the positivity of 'I'xo,^ contributes to the increase of the kinetic energy 
around the point x at scale R. 

Since the flux consists of both the kinetic and the pressure parts, a natural 
question is whether there is a transfer of the kinetic energy from larger scales into 
f?(xo,i?), or perhaps the increase is mainly due to the change in pressure. In gen- 
eral, it is possible that the increase of the kinetic energy around xq is due solely to 
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the pressure part; a simple example being u = (c, c, c) t, p = cx\ + cx 2 + cx$ . How- 
ever, in physical situations where the kinetic energy on the (global) spatial domain 
ft is non-increasing, e.g., a bounded domain with no-slip boundary conditions, or 
the whole space with either decay at infinity or periodic boundary conditions (here, 
we are concerned with the case of decaying turbulence, setting the driving force to 
zero), the increase of the kinetic energy in B(x , R) - and consequently, the positiv- 
ity of $x ,-R ~ implies local transfer of the kinetic energy from larger scales simply 
because the local kinetic energy is increasing while the global kinetic energy is non- 
increasing resulting in decrease of the kinetic energy in the complement. This is 
also consistent with the fact that in the aforementioned scenarios one can project 
the NSE - in an appropriate functional space - to the subspace of divergence-free 
functions effectively eliminating the pressure and revealing that the local flux <J> Xo ,,r 
is indeed driven by transport/inertial effects rather than the change in the pressure. 
(The u = (c, c, c)t, p — cx\ + cx 2 + cx 3 example pertains to a completely opposite 
situation, the kinetic energy is simply uniformly growing over the whole spatial 
domain.) 

Henceforth, following the discussion in the preceding paragraphs - in the set- 
ting of decaying turbulence (zero driving force, non-increasing global energy) - the 
positivity and the negativity of $ Xo ,-R will be interpreted as a transfer of (kinetic) 
energy around the point Xo at scale R toward smaller scales and a transfer of 
(kinetic) energy around the point x at scale R toward larger scales, respectively. 

Completely analogous definitions hold for shells of radii R and 2R. 

We also consider finite time averages of each of the aforementioned quantities. 

Our goal is to obtain a manifestation of the (kinetic) energy cascade in physical 
space, i.e., formulate a condition on B(0,Rq) that would imply that the time- 
averaged energy transfers/cascades to smaller scales across a range of scales (the 
existence of the inertial range). 

A key point here is that we do not assume any homogeneity of the flow; hence 
one can not expect to show that the local fluxes are positive for each individual 
ball B(x, R). The best one can hope for is to prove the positivity of the flux over 
some spatial average. 

We choose to work with a very straightforward spatial average: the arithmetic 
mean of the local fluxes - time-averaged, per unit mass - computed over a family 
of coverings of B(0, R ), the so-called optimal coverings. 

Let K\ and K 2 be two positive integers. A covering {B(xj, R)} ™ =1 of B(0,R n ) 
is an optimal covering (with parameters K\ and K 2 ) if 



and any point x in B(Q,Ro) is covered by at most K 2 balls B(xi,2R). (Optimal 
coverings exist for all large enough K\ and K 2 , the critical values depending only 
on dimension of the space. In M. 3 , we can take K\ = K 2 = 8.) 

Let / be a sign- varying quantity (e.g., the flux density — [(u- V)u + Vp] -u), and 
consider the arithmetic mean of the quantity locally averaged over the (optimal) 
covering elements B(xi,R), 
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(for some < p < 1). 

A revealing observation is that Fr ~ const (R) for all optimal coverings at 
scale R (K\ and fixed) indicates there are no significant fluctuations of sign of 
/ at scales comparable or greater than R. In other words, if there are significant 
fluctuations of sign of / at scale R* , Fr will run over a wide range of values while 
the average is being run over all permissible optimal coverings (determined by K\ 
and K%), for any R comparable or less than R*. 

When there is no change of sign at all, i.e., in the case of a signed quantity (e.g., 
the energy density / = ^|u| 2 , p = 25— I or the enstrophy density / = |V<8>u| 2 , p = 
1), one would then expect that for any scale R, < R < Ro, the averages Fr are 
all comparable to each other. This is in fact true (an easy proof). 

Utilizing the NSE via the local energy inequality - in the mathematical setting 
of suitable weak solutions [17[ [1] - we establish the positivity and near-constancy 
(comparable to vE, where v is the viscosity and E is the average enstrophy over 
B(0,2Rq) x (0, 2T)) of the averaged flux across a range of scales under a very 
simple and natural (in the sense of turbulence phenomenology) condition; namely, 
that Taylor micro-scale ro associated with B(0,Ro) is smaller than the integral 
spatial scale Ro (cf. (|4.13jl ). The larger the gap, the deeper the inertial range. This 
condition is reminiscent of the Poincare inequality on a domain of the corresponding 
size (see Remark |4.2[) ; moreover, the condition in hand would be easy to check in 
physical experiments as the averages involved are very straightforward. In addition, 
the length of the time interval T is consistent with the intrinsic scaling of the model 

( C f. my 

It is interesting to interpret the cascade in the light of the above observation 
regarding the meaning of near-constancy of optimal cover averages. Essentially, for 
any tq < R < Rq (within the inertial range), the flux density does not experience 
significant fluctuations of sign at scale R; the significant fluctuations of sign are 
only possible at the scales substantially smaller than ro, i.e., inside the dissipation 
range. 

The second part of the paper concerns locality of the flux. It is believed (see [16]) 
that the energy flux inside the inertial range of turbulent flows depends strongly 
on the flow in nearby scales, its dependence on the lower and much higher scales 
being weak. The theoretical proof of this conjecture remained elusive. The first 
quantitative results on fluxes were obtained by early 70 's (see [E]). Much later, 
the authors in 15 used the NSE in the Fourier setting to explore locality of scale 
interactions for statistical averages, while the investigation in [5] revealed the lo- 
cality of filtered energy flux under an assumption that solutions to the vanishing 
viscosity Euler's equations saturate a defining inequality of a suitable Besov space, 
i.e., under a (weak) scaling assumption. A more recent work [2] provided a proof of 
the locality of the energy flux in the setting of the Littlewood-Paley decomposition. 

In the last section we prove the locality of the energy cascade - in decaying 
turbulence - in the physical space throughout the inertial range established in 
Theorem 4.1. In particular, considering dyadic shells at the scales 2 k R (k an integer) 
in the physical space, we show that both ultraviolet and infrared locality propagate 
exponentially in the shell number k. 

To the best of our knowledge, the condition (|4.13j) is presently the only condition 
(in any solution setting) implying both the existence of the inertial range and the 
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locality of the energy flux. Moreover, it does not involve any additional regularity 
or homogeneity /scaling assumptions on the solutions to the NSE. 

2. PRELIMINARIES 

We consider three dimensional incompressible Navier-Stokes equations (NSE) 

-|u(t, x) - vAu(t, x) + (u(t, x) • V)u(t, x) + Vp(t, x) = 0, 

V • u(i,x) = 

where the space variable x is in R 3 and the time variable t is in (0, oo). The vector- 
valued function u and the scalar-valued function p represent the fluid velocity and 
the pressure, respectively, while the constant v is the viscosity of the fluid. 

Since our goal is to investigate local fluxes in the physical space, the class of suit- 
able weak solutions (see [HE]) will provide an appropriate mathematical frame- 
work. 

Definition 2.1. Let SI be an open connected set in R 3 . We say that (u,p) is a 
suitable weak solution on (0, oo) x SI if 

(a) u€ L oo ((0,oo),L 2 (Sl) 3 )nL 2 ((0,oo),H 1 (Sl) 3 ) and p e L 3/2 ((0, oo) x SI); 

(b) the NSE i2. 1)) are satisfied in the weak (distributional) sense; 

(c) the local energy inequality is satisfied: for any <f> € X>((0, oo) x SI), <p > 
we have 

(2.2) 

2v J J |V ® u\ 2 (j)dxdt < J J \u\ 2 (d t (t) + vA(t>)dxdt + J J (\u\ 2 + 2p)u ■ \7 (f) dx dt 

where 2?((0,oo) x 0) denotes the space of infinitely differentiable functions with 
compact support in (0, oo) x St. 

The existence of the suitable weak solutions in the case where SI = M 3 and 
the external force is zero, given a divergence- free initial condition in i 2 , was first 
established in [17] . See also [1] [14] for more general results related to existence and 
regularity properties of the suitable weak solutions. 

A solution to the NSE on (0, oo) x SI is called regular if its H 1 norm is bounded on 
(0, T) for any T positive. Given appropriate boundary conditions, this implies that 
the solution is infinitely differentiable (in fact, analytic) in both space and time and 
so it is a classical physical solution. In particular, the local energy equality holds 
( (|2.2|) becomes an equality). The smoothness of the suitable weak solutions to the 
NSE is still an open problem, and the best result in this direction reads that the 
one-dimensional (parabolic) Hausdorff measure of the singular set in (0, T) x SI is 
zero pQ (outside the singular set, a suitable weak solution is infinitely differentiable 
in the spatial variables). 

In what follows, we consider 

(2.3) R Q > such that B(0, 3i? ) C SI 

where -6(0, 3Rq) denotes the ball in R 3 centered at the origin and with the radius 

Let 1/2 < 5 < 1. Choose tp € T>(B(0,2Rq)) satisfying 

| V^o| Co [Aj^pj Co 



(2.4) <*,<!, V^lon5(0,i?„), ™ < ^<^> 
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region containing the rest of supp(^) 




FIGURE 1. Regions of supp(V>) in the case B(xq,R) <£ 5(0, i? ), 
cross-section. 



For a T > (to be chosen later), x G B(0,R ) and < R < R define 
4> = 0x o ,T,_R(i, x) = r)(t)ip(x) to be used in the local energy inequality (|2.2j) where 
^ = 'Tr(^) arl d V' — ^xo,-r( x ) are refined cut-off functions satisfying the following 
conditions, 

(2.5) Tje2?(0,2T), < 77 < 1, V = 1 on (T/4, 5T/4), ^J<^; 

if B(x ,R) C 5(0, i? ), then V> G V(B(x Q ,2R)) with 
(2.6) 

|W| / Co |At/;| „ Co 



0<V><V> , V> = 1 on B(x ,R)n B(0,R Q ), 



< 



< 



ip s ~ R 1 iP 26 - 1 " i? 2 
and if B{x ,R) £ B(0,R ), then -0 G V(B(0,2R a )) with 0> = 1 on B(x ,R) D 
B(0,Ro) satisfying, in addition to (|2.6p . the following: 

■0 = "00 on the part of the cone in R 3 centered at zero and passing through 



(2.7) 



and 



5(0, Rq) n B(xq,R) between 5(0, R ) and 5(0, 2R ) 



ij) = on B(0, Rq) \ -B(xrj, 2R) and outside the part of the cone in M 3 

(2.8) centered at zero and passing through 5(0, Rq) n £>(xo, 2i?) 

between 5(0, Rq) and 5(0,2i? )- 

Figure [1] illustrates the definition of ip in the case B(xo,R) is not entirely con- 
tained in B(0,R ). 



Remark 2.1. The additional conditions on the boundary elements (|2.7|) and (|2.8j) 
are necessary to obtain the lower bound on the fluxes in terms of the same version 
of the localized enstrophy E in Theorems 14.11 and 15.21 (see Remarks 14.41 and 15. 5|) . 
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3. Localized Energy, Enstrophy and Flux; Ensemble Averages 

Let xo £ B(0,Rq) and < R < Rq. Define localized energy, e, and enstrophy, 
E, at time t - all per unit of mass - associated with i?(xo, R) by 

(3-1) e XOifi (i)= [ huf^*- 1 dx , 



2 

(3.2) E X0 , fl (t) = j |V®u| 2 0dx 

(for some | < 8 < 1). 

The total - (kinetic) energy plus pressure - flux through sphere S(xq, R) is given 

by 

J (i |u| 2 +p)u ■ nds = J [(u ■ V)u + Vp] • urfx 

where n is an outward normal. Considering the NSE localized to B(xq, R) leads to 
a localized version of the flux, 



(3.3) 



$x 0> flW = / (\\ U \ 2 +P) U - V0dx= -J [(u- V)u + Vp] -u0dx 



where <p — with n and "0 as hi (I2.5H2.6|) . Since ip can be constructed such that 
= is oriented along the radial directions of B(xq, R) towards the center of 
the ball xo, <l>(xo,i?) represents the flux into i?(xo,i?) through the layer between 
the spheres S(xq, 2R) and 5(xo, R) (in the case of the boundary elements satisfying 
the additional hypotheses (|2.7[) and (|2.8|) . ip is almost radial and the gradient still 
points inward). 

For a quantity 6 Xi _r(^), t € [0, 2T] and a covering {-B(xj, i?)}i=i, n of B(0,R ) 
define a time-space ensemble average 

(3-4) <e>fl = i /^E^Q*i,fl(*)^- 

z— 1 

Denote by 



(3.5) 


ej? = 


(e x ,fl(*)>fl 


(3.6) 


E R = 


(Ex,B(t))B 


(3.7) 


$il = 





the averaged localized energy, enstrophy and inward-directed flux over balls of ra- 
dius i? covering B(0, Rq). 

Also, introduce the time-space average of the localized energy on B(0, Ro), 

(3-8) e = i / 4 e o,%(*) * = 44* // ilul^iT 1 rfxdt 



T J Rl ' T RlJJ 2' 

and the time-space average of the localized enstrophy on £?(0,i?o), 

(3.9) E = ^/ Rj E ^ Ro{t)dt = l ; R^ JJ l V ® u | 2 ^ dxd< 
where 

(3.10) 0o(t,x)=»7(t)^o(x) 
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with ip defined in (|2T4|) . 

Finally, define Taylor length scale associated with B(0,Rq) by 

/e \ V2 

(3-11) to = (I) . 

Note that the possible lack of regularity may produce additional loss of energy, 
resulting in anomalous energy dissipation and the loss of flux leading to the strict 
inequality in (|2.2p . Let us mention here that in the turbulence literature the term 
'anomalous dissipation' is usually utilized in the context of the possible energy dissi- 
pation due to the (possible) singularities in the 3D Euler equations (the observation 
originally made by Onsager); for rigorous results on Onsager's conjecture on the 
energy conservation in the Euler equations see, e.g., [3], and a recent work [2J. 

Denote by R the loss of flux due to possible singularities in [0, 2T] xB(xo, 2R) , 

r (i|u| 2 +p)u-V^xdi 

® u\ 2 <j)dxdt - i ff \u\ 2 (d t (j) + uA(/))dxdt 



C-S.12) " 2 

= // / / IV n I J , w/x ill - ^ 

where <j) = r/ip with n and ip as in (12.5[) and (l2.6H2.8p . In particular, denote by 
(J)^ = the loss of flux due to singularities in [0, 2T] x 5(0, 2R ). 

We will also consider the time-space ensemble averages of these anomalous fluxes, 



1 " 1 1 

(3.13) $^ = -^-— 



TR 3 ~ Xi ' R 



i=i 

Note that due to (I2.2[) . all the anomalous fluxes are nonnegative, 
(3.14) $£U>0, $oo>0, $^>0; 

they are all zero provided the equality holds in (I2.2j) inside [0, 2T] x B(0, 2Ro). In 
particular, the anomalous fluxes are all zero provided the solution in view is regular 
on [0,2T] x B(0,2i? ). 

Consequently, the total localized flux into i?(xo, R) over interval [0, 2T] , including 
the (loss of) flux due to the possible loss of regularity, is 

(3-15) * X0 ,« = J ®x , R (t)dt-$Z,M 

and the time-space ensemble average of this flux at scales R over time T is 

(3.16) *^ = $ fl -$£. 

We will refer to ^xo.-R an d as the modified flux over [0, 2T] into B(xo,i?) and 
the (time-space ensemble) averaged modified flux at the scale R, respectively. 

Let Ki,K2 > 1 be two positive integers (independent of R,Ro, and any of the 
parameters of the NSE). 

Definition 3.1. We say that a covering of 13(0, Rq) by n balls of radius R is 
optimal if 



(3.18) 



any x G -B(0, Rq) is covered by at most K2 balls i?(x;, 2R) 
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Note that optimal coverings exist for any < R < Ro provided K\ and K2 are 
large enough. In fact, the choice of K\ and K2 depends only on dimension of the 
space; in R 3 we can choose K\ = K2 = 8. 

Henceforth, we assume that the averages are taken with respect to optimal 
coverings. 



Lemma 3.1. If the covering {i?(x.;, J2)}i=i, n of B(0, Rq) is optimal then 

1 1 

where K > is a constant depending only on K 2 and dimension of the space R 3 . 



(3-19) 



Proof. Let {x^ } be a subset of {x.i}i=i, n such that interiors of the balls B(x^ , 2R) 
are pairwise disjoint. Using (|3.12l) . we obtain 



(3.20) 

and 
(3.21) 



(-|u| 2 +p)u- Vtpodxdt - $ c 



= v J J |V ® u\ 2 (f)odxdt - J J \u\ 2 (d t (j)o + vAfio) dxdt 

i|u| 2 + P )u.v(^^.)dxdi -j>~ m 

3 3 

= vjj |V®u| 2 g>,)d*it-l If |u| 2 [S t (^ ^J + i/A^ 

3 3 3 

where (f>o = rjipo and </>j. = rpp^ with rj as in (|2.5p . i/)q &s in (|2.4p and ij)^ a test 
function corresponding to £?(xj. ,i£) satisfying (j2.6H2.8p . 
Note that the definitions of 4>o and 4>i j imply 

4> = 0o - hi ^ > 

3 

hence, by the local energy inequality 
(3.22) 



(i |u| 2 +p)u ■ \70dxdt 



~ V IS ' V<8) U \ 2 ^ dxdt ~\ II \ u \ 2 ( d t<i> + vA<t>)dxdt . 
If we add relations (|3.2ip and (|3.22p and then subtract (|3 . 20[) we obtain 
(3-23) $00 

3 

Let £ be a cubic lattice inside B(0,Rq) with the points situated at the vertices 
of cubes of side R/2 (Note that this lattice can be chosen such that the number of 
points in it is between 2 3 (i? /i?) 3 and (47r/3)2 3 (i? / J R) 3 ). 

Since the covering {_B(xi, R)} is optimal, each point in C is contained in at most 
K2 balls. Moreover, any ball in the covering will contain at least one point from 
the lattice. 

If CJ is sub-lattice of C with points at vertices of cubes of side 4R, then the 
interiors of balls of radius 2R containing different points of £' are pairwise disjoint, 
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and thus if we denote by -B(x; , R) a ball from the covering {B(xi, R)} containing 
the point p e £, by 

Note that for each point p E C there are at most if 2 choices for i?(xj , i?). So 

if 2 $oo > X! • 

i:B(x»,il)n.C'#0 

Clearly £ can be written as a union of 8 3 = 256 sub-lattices £' k , k = 1, . . . , 256, 
each C' k having the same properties as C . Thus, 

n 
i=l 

Consequently, 

$n=--^-V$~n< 8 3 if 2 --^-$oo < 8 3 if 2 -^r*oo 
fl TR?n^ X " R ~ T R 3 n T R% 

1=1 u 

where the last inequality is due to n satisfying (|3.18[) . □ 

According to the lemma, the time-space ensemble averages taken over the 
optimal coverings at the scale R are bounded, independently of i?, by the average 
loss of flux due to possible singularities inside 5(0, 2i?o). 

4. Energy Cascade 
Let {.B(xj, i?)}i=i,7j be an optimal covering of -B(0,i? ). 

Note that the local energy equality ([3.1 2|) and the definitions of E#, <J>/j and <£>^ 
( ([311), ((3l7j) and flam) ) imply 

(4.1) y H = $ R -$<g= v E R -±-f2±;-L f[±\u\ 2 (d t <f> i + vAcf> i )dxdt 

i—l 

where <pi — rjipi and ipi is the spatial cut-off on B(xi,2i?) satisfying ([2.5H2.8p . 
If 

R 2 

(4.2) T> " 



1/ 



then for any < R < R , 



1 *"/, ^ ,P Q J.2S-1 



U l )t\ = \m^\<C a -r 1 A ^<v^. 



hence. 



C 

> vE R - v— e R . 



Co i.25-1. 



The optimality conditions (j3.17[) and (13.181) paired with (j2.6H2.8p imply 

(4-3) E fl > -Le 

ivi 
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and 
(4.4) 

Consequently, 

(4.5) 



leading to the following proposition. 
Proposition 4.1. 

(4.6) *r>c x ve(i-c^ 

with C\ = 1/ K\ and c 2 = CqK\K 2 (provided conditions \3. 1 TfOTTty) are satisfied). 
Suppose that 



(4.7) 



C 2 



1 /2 

for some < 7 < 1. Then, for any R, (c 2 /7J t < R < Rq, 
(4.8) >ci(l- 7 2 )^E = c , 7 ^E 

where 



(4.9) 



Co.- 



ci(l-7" 



1-7 2 
#1 



To obtain an upper bound on the averaged modified flux, note that for optimal 
coverings, in addition to (|4.4p . 



(4.10) 

Hence, (|4.1|) implies 



E R < K 2 E 



C 1 
*n < vE H + v Jfi e R <vK 2 E + i>C K 2 — e. 

If the condition (|4.7|) holds for some < 7 < 1, then it follows that for any R, 
(c2 1/2 /i)t <R<Ro, 



(4.11) 
where 
(4.12) 



,Tr I ' C C0-^27 2 c ^ c 

Wj{ < ^A2E + f E < C\ a Vt. 



Cl, 7 = #2 



1 



g 7 2 

C2 



f'2 



= #2 



1 



7 



^1^2 



Thus we have proved the following. 
Theorem 4.1. Assume that for some < 7 < 1 

(4.13) ro<c 7 i?o, 
where 

(4.14) c= — 1 = . 
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7 2 



Then, for all R, 

(4.15) — t <R<Rq, 

the averaged modified flux "P/? satisfies 

(4.16) c , 7 ^E <^ R < Cx^vE 
where 

1 -7 2 r 

(4-17) c , 7 = — i- , c 1)7 = Jf 2 1 + , . 

and i/ie average (•)# is computed over a time interval [0, 2T] wii/i T > R\jv and 
determined by an optimal covering of B(Q,Rq) (i.e., a covering satisfying Jff._?7| ) 

and fug) ;. 

Remark 4.1. As noted in the introduction - in the case the global energy is non- 
increasing - the theorem provides a sufficient condition for the energy cascade; i.e., 
a nearly constant nonlinear transfer of time-averaged (kinetic) energy to smaller 
scales across the inertial range defined by (|4.15p . 

More precisely, since we are working with weak solutions, the expression for the 
rate of change of local kinetic energy (|1.1[) morphs into 

(4.18) — J J i|u| 2 t = ^xo.if + viscous terms. 

The interpretation remains the same, the only differences being that the time- 
derivative of the local kinetic energy is taken in the sense of distributions and the 
flux got replaced with the modified flux to account for possible singularities. 

Note that the averages are taken over finite-time intervals [0, 2T] with T > Rq/v 
(see (|4.2[) ). This lower bound on the length of the time interval T is consistent 
with the picture of decaying turbulence; namely, small v corresponds to the well- 
developed turbulence which then persists for a longer time and it makes sense to 
average over longer time-intervals. 

Remark 4.2. In the language of turbulence, the condition (|4. 13[) simply reads 
that the Taylor micro-scale computed over the domain in view is smaller than the 
integral scale (diameter of the domain). 
On the other hand, (|4.13p is equivalent to 

1 f /Y.I2-M-1 , 7 2 n2 1 



j u \u\'4tf- 1 dxdt < ^ Si- JJ |V®u|'0odxdt 

which can be read as a requirement that the time average of a Poincare-like in- 
equality on -8(0, 2i?o) is not saturating; this will hold for a variety of flows in the 
regions of active fluid (large gradients). 

Remark 4.3. Since $^ > 0, all the lower bounds on the averaged modified fluxes 
hold for the usual averaged fluxes in particular, 

(4.19) $fl>c^E (l-c 2 ^ 



R 2 

and consequently, provided (I4.13|) and (|4.15[) hold, 
(4.20) $ fl > cq j7 i/E . 
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Also, if a solution u is such that the (|2.2|) holds with equality and in particular, 
if the solution is regular, then all the estimates, including (|4.16p . hold for the usual 
averaged flux 

Remark 4.4. If we do not impose the additional assumptions (|2.7[) and (|2.8[) for 
the test functions on the balls B(xi,R) <£_ B(0,Rq), then the lower bounds for 
in (|4.6|) and (|4.16j) will hold with E replaced by the time-space average of the 
non-localized in space enstrophy on B(0, Rq), 

it 

E ' = k\h J i v ® u iV xdi - 

°B(0,flo) 

This is the case because the estimate (14.3[) gets replaced with 

E R > J-E' . 
~ K\ 

5. Locality of the averaged flux 

Let Xo G B(0, Rq), < R2 < R\ < Ro- In order to study the flux through the 
shell A(x , Rx, R2) between the spheres S(x ,R 2 ) and S(x ,Rx), in what follows, 
we will consider the modified test functions (f> — 4> Xo ^t,r 1 ,r 2 (£> x ) — 7 7(OV'( x ) to 
be used in the local energy inequality (|2.2|) where 77 = ?7t(£) as in f|2.5f) and ijj £ 
V(A(x ,2Rx,R 2 /2)) satisfying 

(5.1) 0<V.<^o, ^ = lon^(x 0)J Ri, J R 2 ) ) ^ < < § 
where ipo is defined in (|2.4p and 

(5.2) R = R{Rx,R 2 ) = min{i? 2 ,i?i - R 2 } ■ 

Use the <fi above to define the localized time-averaged flux through the shell 
between the spheres <S'(xo,i?2) and S(jx.q,Rx) as 

(5-3) ^.jfe.jfe = - /7(^|u| 2 +p)u-V<Hxdi ■ 

If Ri R2 is the anomalous flux inside A(x , 2i?i, R 2 /2), i.e., if $x ,fli,fl 2 sat ~ 
isfies 



(5.4) 



(i|u| 2 + P )u-V0dxdt -*~ KliJJa 

|V<E> u| 2 0dxd< - i ^ \u\ 2 (d t ct) + vA(t))diidt 



then define the time average of the modified localized fluxes through the shells 
A(xo,Rx,R 2 ) by 

(5-5) *xo,iii,fla = *x D ,fli,B 2 - 7^®x ,R 1 ,R 2 ■ 

As already mentioned, the modified fluxes can be viewed as total fluxes including 
the (loss of) flux due to possible singularities inside the shell. Also note that the 
local energy inequality (I2.2p implies 

(5-6) > 0. 
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Define the time-averaged energy and enstrophy per unit of mass in the shell 
between the spheres 5*(xo,i?2) and S(yio,Ri) by 

(5-7) 1 

Exo.Bi.Ba = j, / / |V ® u| 2 0dxdt ; 

then, 

1/2 



( 5 - 8 ) Txo,fli,fi 2 



E Xo ,fl:i,i?2 



is the Zoca/ Taylor length scale associated with the shell A(xo, Ri, R2). 
Note that 

(5.9) u\A4>\ = v\r,M\ < Av^l < "S^ -1 



and 



1 S„i, ^ ,25-1 



(5.10) W = M|<C -^<^ 

provided 

(5.1!) 

Hence, (|5.4I) implies that for any xo € B(0,Rq) and any < R2 < Ri < R, 

Co 

*x 0) fli,B 2 > ^E XOiiil)fi2 - v-~^e XQ:Rl . R . 2 
(5-12) / 

r xo,Bi,fl 2 



leading to the following proposition. 

Proposition 5.1. Let < 7 < 1. Then, for any shell ^4(xo, R\, R2) satisfying 

(5.13) t XOi r 1: r 2 < -^j^R 

Co 



with R defined by 15. ty . 

(5-14) ^xo.Bi.ifa > vE^r u r 2 {1 - 7 2 ) . 

Similarly, utilizing f|5 .4[) again, we obtain an upper bound 

n Co 
Vx ,fli,B 2 < fEx ,fli,fl 2 + ^7v7 e x ,Bi,B 2 ■ 

If the condition (15.131) holds for some < 7 < 1, then it follows that 
(5.15) ^x ,Bi,R 2 < ^E X0)Rl)fl2 + 7 2 ^E X0i _ Rl ^ 2 < (1 + 7 2 ) vE^ 0tRuR2 ; 

thus, we have arrived at our first locality result. 
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(5-16) T XOtRuRa < —y^R 



Theorem 5.1. Let < 7 < 1, x € B(0, R ) and < R 2 < Ri < R . If 

j_ 

'0 

i? defined by t5.ty) . then 

(5-17) (1 - 7 2 ) ^E X0)Rl)i?2 < * X0)fll ,H 2 < (1 + 7 2 ) ^E X0 , fllifl2 

where the time average is taken over an interval of time [0, 2T] with T > R^/v. 

Remark 5.1. The theorem states that if the local Taylor scale t^^^, asso- 
ciated with a shell A(xo, Ri, R2), is smaller than the thickness of the shell R (a 
local integral scale) , then the time average of the modified flux through that shell 
towards its center Xo is comparable to the time average of the localized enstrophy 
in the shell, E Xo ,r 1: r 2 . Thus, under the assumption (|5.16|) the flux through the shell 
^4(xo, Ri, R2) depends essentially only on the enstrophy contained in the neighbor- 
hood of the shell, regardless of what happens at the other scales, making (|5. 13|) a 
sufficient condition for the locality of the flux through ^4(xo, Ri, R 2 )- 

Remark 5.2. Similarly as in the case of condition (|4.13|) . we can observe that 
condition (|5.16|) can be viewed as a requirement that the time average of a Poincare- 
like inequality on the shell is not saturating making it plausible in the case of intense 
fluid activity in a neighborhood of the shell. 

Remark 5.3. Since <&°° „ B > 0, all the lower bounds on the modified fluxes 
hold for the usual fluxes c E> X o.i?.o,.Ri ! m particular, we have 

(5.18) ** ,Rr,R 2 > ^E X0 ,fl 1)fl2 (\ - Co T *oA,^ 

and, provided (|5.16j) holds, 

(5-19) *x ,iio,fli > (1 - 7 2 ) ^E X0>Kliil2 . 



Also, if a solution u is such that the (|2.2j) holds with equality and in particular, 
if the solution is regular, then all the estimates, including (|5.17j) . hold for the usual 
averaged flux ^xo,^,^. 

In order to further study the locality of the flux, we will estimate the ensemble 
averages of the fluxes through the shells A(xi,2R, R) of thickness R = R. Since 
we are interested in the shells inside B(0,Rq), we require the lattice points to 
satisfy 

(5.20) B( Xi ,R) c B(0,Rq) . 

To each A(x,;, 2R, R) we associate a test function <fii = rjtpi where r\ satisfies ()2.5[) 
and ipi satisfies the following. 

If A(xj, 2R, R) C B(0, R Q ), then ip % e X>(A(x 4 , 4R, R/2)) with 
(5.21) 

< i/>i < Vo, 4>i = 1 on A^, 2R, R) f] B(0, R Q ), < ^, -^g < ^ , 

and if A(Xi,2R,R) (£ 5(0, R n ) (i.e. we have B(x t ,R) C B(0,R ) and B(x u 2R) \ 
B(0,R a ) ^ 0), then ^ £ V(B(0,2R )) with ^ = 1 on A(x ,2i?,i?) n B(0,R a ) 
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4>i = </>o 



region containing the rest of supp(j^j) 




Figure 2. A cross-section of regions of supp(-0i) in the case 
A(x h 2R,R) t B(0,Ro). 



satisfying, in addition to (I5.2ip . the following: 

■0i — *Po on the part of the cone in R 3 centered at zero and passing 



(5.22) 



and 



through 5(0, R ) H B(x i; 2R) between 5(0, R ) and 5(0, 2R Q ) 



ipi — on B(0, i?o) \ A(x.i, 4R, R/2) and outside the part of the 

(5.23) cone in R 3 centered at zero and passing through 5(0, i?o) H B(xi,4R) 
between 5(0, _R ) and S(0,2R ). 

Figure [5] illustrates the definition of if>i in the case A(x.i,2R, R) is not entirely 
contained in B(0 7 R ). 

Similarly as in the previous section, we consider optimal coverings of B(0, Rq) 
by shells {A(x.i, 2R, R)}i = i. n such that (|5.20|) is satisfied, 

(5.24) ^)%„< A ,(^ 3 
and 

(5.25) any x G B(0, R ) is covered by at most K 2 shells A(xi, 4i?, R/2) 
Introduce 

1 " 

hR,R = - V" e x 



(5.26) 
(5.27) 



,2R,R 



i = l 



1 - 

£-2R,R = — / E : 



Xi,2R,R. 
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and 

1 - 

(5-28) $ 2 R,R = - V <S>x z ,2R,R , 

i=l 

the ensemble averages of the time-averaged energy, enstrophy, and flux on the shells 
of thickness R corresponding to the covering {A(x.i, 2R, R)}i=i, n ■ 

The ensemble average of the time-averaged modified flux on shells of thickness 
R is then defined by 

1 ™ 

(5.29) ^ 2 R,R = - **i,2fl,fl = $>2R,R ~ ®Tr,R 



n 
i=i 



where 



(5-30) ^Tr,r = ^J2 ^Z2R,r 

i=\ 

is the ensemble average of the time-averaged anomalous fluxes through the shells 
of thickness R inside B(0,Rq). 

An argument analogous to the one in Lemma 13.11 implies that if the covering of 
B(0, Rq) is optimal, then 

(5-31) $? RtR < A^Soo ■ 

Taking the time ensemble averages in (|5.4|) and applying the bounds (|5.9ti5. 10]) . 

we arrive at 

Co 

(5.32) ^2R,R > V^2R,R ~ V— e 2 R,R ■ 

If the covering is optimal, i.e., if (|5.20|) and (|5.24H5.25[) hold, then 

1- 1 ( R ^ 



(5-33) E 2 r,r > -E > — — 

n Ki \R 

and 

(5.34) ~e 2RiR <^l~e<K 2 

n \R , 

where 

(5.35) E = ^ J J |V ® u| 2 o dxdt — RqE 

is the time average of the localized enstrophy on B(0, Rq) and 
(5-36) S= 5?// \ u \ 2 ^o 5 ~ ldxdt = R o^ 

is the time average of the localized energy on B(0, i? ) with cj) Q is defined by (|3.10p . 
Let us note that 

1/9 / ~ \ 1/2 

/ e \ 12 I e \ 

(5.37) „_y ^ 

Utilizing (jOSf . ([5734]) and ([537]) in the inequality (f5T32|) gives 

(5.38) ^^(D^K 1 "^ 2 ^ 
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implying the following result. 

Proposition 5.2. Assume that the condition ^4.13ty holds for some < 7 < 1. 
Then, for any R satisfying l[4-15\ ), 

R^ 3 



(5.39) c , 7 l-j < * 2i? ^ 

ZioZds im'i/i c and co l7 defined in \4-lJ$ and J^.i7| ). 

Taking the time ensemble averages in the localized energy equality (|5.4[) again, 
this time looking for an upper bound, yields 

t ^ Co ~ 

^2fl,fl < Vt 2 R,R + V — e 2 R,R ■ 

If the covering {A(x i} 2R, R)}i=i, n of i?(0, R ) is optimal then, in addition to (|5.34|l . 

(5.40) E 2R;R < — E < K 2 (-f 

n \Rn 



hence, 



r 1 R \ T , Cn / R X 



7?o y i? 2 y i?Q 

As long as R is inside the inertial range delineated in (|4.15l) . 

e < c 2 7 2 E 2 E, 

leading to 

R\ ~ Tr „ / \ 9 n~ / R ^ 



* 2 r,r<»K 2 1-\ E + vK 2 C Q \—j ^fE = c li7 l-l i/E. 

Collecting the bounds on ^2/?,^ we establish the following. 

Theorem 5.2. Assume that the condition fr4-13\ ) holds for some < 7 < 1. Then, 
for any R satisfying the ensemble average of the time-averaged modified flux 

through the shells of thickness R, ^ 2 R y R, satisfies 

(5.41) c Qn (Jj^J vE < V 2R . R < ci, 7 (-^j vE 

where c, Co, 7 , and ci i7 are defined in ^4-14\ ) and j^.l7[ ) and the average is computed 
over a time interval [0, 2T] wrai/i T > R^jv and determined by an optimal covering 
{A(xi,2R,R)} i=hn ofB(0,Ro) (i.e. satisfying flTMj) . JO^ , and {Q^ J. 

Note that if 

^2R,R = J^ i &2R,R 

denotes the ensemble average of the time- space averaged modified flux through the 
shells of thickness R then, dividing (|5.41[) by R 3 , we obtain the following. 

Corollary 5.1. Under the conditions of the previous theorem, 

(5.42) c , 7 ^E < ^ 2RtR < c hl isE . 
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Theorem 15.21 allows us to show locality of the time-averaged modified flux under 
the assumption (|4.13[) . Indeed, the ensemble average of the time-averaged flux 
through the spheres of radius R satisfying (14. 15)) is 



According to Theorem 14.1 



R\ 3 ~ ~ ( R^ ' 



On the other hand, the time ensemble average of the flux through the shells between 
spheres of radii R2 and 2i? 2 , according to Theorem 15.21 is 

Z? \ / E? ^ 



C0, 7 ( J ^2R 2 ,R 2 < Cl, 7 I ) f E 



c 0, 7 / ^2R 2 ,R 2 Cl, 7 / R2 



Consequently, 

(5-43) ^(77) - 

Cl, 7 \ R J C , 7 V R 

Thus, under the assumption (|4.13p . throughout the inertial range given by (|4.15p . 
the contribution of the shells at scales comparable to R is comparable to the total 
flux at scales R, the contribution of the the shells at scales R2 much smaller than R 
becomes negligible (ultraviolet locality) and the flux through the shells at scales R2 
much bigger than R becomes substantially bigger and thus essentially uncorrelated 
to the flux at scales R (infrared locality). 

Moreover, if we choose R2 = 2 k R with k an integer, the relation (|5.43[) becomes 

(5.44) ^l2 3fc < *2 fc+ _ 1J *- 2fc * < 2hl 2 3k 

Cl, 7 fy R Co, 7 

which implies that the aforementioned manifestations of locality propagate expo- 
nentially in the shell number k. 

In contrast to (|5.43[) . since e = i?j]e, E = RqE, ^2R 2 ,r 2 = R2^2R 2 ,r 2 and 
V R = R 3 V R , 

(5.45) ^2 < < £l£L , 

Cl, 7 C ,7 

i.e., the ensemble averages of the time- space averaged modified fluxes of the flows 
satisfying (|4. 13[) are comparable throughout the scales involved in the inertial range 
(I4.15|) which is consistent with the universality of the cascade. 

We conclude this section by noticing that the remarks similar to those at the 
end of section 2] can be applied here. Namely we have the following. 

Remark 5.4. Since R > 0, then all the lower bounds on modified fluxes hold 
for the usual fluxes <&2R,r and <&2R,r] in particular, we have 

R\ 3 =. (. T?' 



(5.46) $2R, R > ci f-j v£ (l - c 2 J, 
and, provided (|4~T3l> and ([47T5]) hold, 

(5.47) $ 2 R,R > c ,7 (J^j i/E. 
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Also, if a solution u is such that (|2.2p holds with equality and in particular, if 
the solution is regular, then all the estimates, including (|5.41|) and (|5.42|) . hold for 
the usual averaged fluxes &2R,r and &2R,r = &2R,r/R 3 - 

Remark 5.5. If the additional assumptions ()5.22|) and (|5 . 23f) for the test functions 
on the shells A(x,, 2R, R) which are not contained entirely in B(0,Ro) are not 
imposed, then the lower bounds in (|5.38p and (15.41)) hold with E replaced by the 
time average of the non-localized in space enstrophy on B(0, Ro), 

2T 

E' = U J \V ®u\ 2 T]dxdt = RlE' . 

B(0,Ro) 

This is the case because the estimate (|5.33[) gets replaced with 

E 2RR >J-(^-) 3 E>. 
2RR ~ Ki\RoJ 

Also, the estimates (|5.43|) and (|5.45j) will contain the terms E'/E(= E'/E) in the 
lower and E/ E' in the upper bounds. 
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